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ABSTRACT 

Large differences between the properties of the known sample of cataclysmic variable 
stars (CVs) and the predictions of the theory of binary star evolution have long been 
recognised. However, because all existing CV samples suffer from strong selection 
effects, observational biases must be taken into account before it is possible to tell 
whether there is an inconsistency. In order to address this problem, we have modelled 
the impact of selection effects on observed CV samples using a Monte Carlo approach. 
By simulating the selection criteria of the Palomar-Green (PC) Survey, we show that 
selection effects cannot reconcile the predictions of standard CV evolution theory with 
the observed sample. More generally, we illustrate the effect of the biases that are 
introduced by magnitude limits, selection cuts in U — B, and restrictions in galactic 
latitude. 

Key words: binaries - stars: dwarf novae - novae, cataclysmic variables. 



1 INTRODUCTION 

Cataclysmic variable stars (CVs) are semi-detached binary 
stars with orbital periods (Porb) typically of the order of 
hours, consisting of a white dwarf primary accreting from a 
companion th at is usu a lly a l ate-type, approximately main- 
sequence star. IWarnerl i| 19951 ) gives a comprehensive review 
of the subject. 

An important question in current CV research concerns 
the secular evolution of these systems. The evolution of CVs 
is of interest in its own right (it determines the size and prop- 
erties of the Galactic CV population), but also because the 
processes thought to drive it are equally relevant to all close 
binary systems containing a late-type component. Binary 
star evolution is driven by loss of orbital angular momentum 
(J) through gravitational quadrupole radiation, and, in at 
least some cases, t he much more efficient mechanism of mag- 
netic braking (e.g. Paczynski 1981; Paczvhski fc Sienkiewiczl 



ll98ll : IVerbunt fc Zwaanlll98ll : lRappaport et al.lll983t ) 
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(e.g. lWarne"rlll995l) . where Mi and M2 are the white dwarf 
and secondary mass respectively, Mi and M2 are their time 
derivatives, and q — M2/M1 is the mass ratio. Equation [T] 
implies that the orbital period distribution of CVs is a useful 
indicator in the study of their evolution — loss of angular 
momentum and the resulting mass transfer changes P or b- In 
practise, P or b is also the only parameter accurately measured 
for a large number of systems. Most studies of the evolution 
of CVs are therefore aimed at accounting for the observed 
Porb distribution. This distribution is shown in the left hand 
panel of Fig. [1] Its most striking features are the period 
gap (a pronounced drop in the number of systems at 2 h < 
Porb i$ 3 h), and the period minimum (a sharp cut-off in the 
period distribution of hydrogen rich CVs at about 80 min). 

The period gap is usually expl ained as resultin g from 
disrupted magnetic braking (e.g. Robinson et al.l Il98ll ; 
iRappaport et al.lll983l ; ISpruit fc Ritterill983l ). In a CV ap- 
proaching the upper edge of the period gap, M2 is approach- 
ing the value where the internal structure of a low-mass 
main-sequence star changes from a deep convective enve- 
lope to fully convective. When the secondary loses its radia- 
tive core, its magnetic field structure is proposed to change, 
causing the magnetic braking mechanism to become less ef- 
ficient, or to cease operating altogether. Angular momentum 
loss then proceeds at the much lower rate allowed by gravi- 
tational radiation. Because of its history of rapid mass loss, 
the secondary is at this point out of thermal equilibrium, 
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Figure 1. The orbital period distribution of non-magnetic, 
hydrogen-rich CVs (left hand panel). Periods are taken from 
I Putter fe Ko lb (2003). The right hand panel shows the same dis- 
tribution divided into two magnitude bins, V > 16.5 (black) 
and V 16.5 (dar k grey). The magnitudes are magi from 
iRitter fe Kolbl J2003T) . and are actually not in all cases V. We 
chose V = 16.5 as a divide because this is roughly the magnitude 
where it becomes difficult to measure P or b from radial velocity 
variations with a 2-m class telescope. The period gap is indicated 
by light grey shading of the range 2.2 h < P or b < 2.8 h in both 
panels. 



which it now regains by shrinking in radius, thereby losing 
contact with its Roche lobe, and only regaining it when P or b 
has decreased to approximately 2 h. 

The disrupted magnetic braking model is broadly con- 
sistent with observational estimates of mass transfer rates 
(M) in CVs — these are on average much higher for long- 
peri od systems (P or b ^ 3 h) than for systems at P or t < 3 h 
(e.g. |Patterson|[l983 ) . That magnetic braking should be dis- 
rupted has, however, been challenged bot h on theoretical 
jTout fc Pringlel Il992l ) and observational l|Andronov et al.l 
120031 ) grounds. Nevertheless, disrupted magnetic braking re- 
mains the conventional explanation for the period gap, and 
a cornerstone of CV evolution theory. 

The period minimum is a consequence of the response 
of a seco ndary of very small mass to continuing mass 
loss (e.g. iPaczvnskil Il98ll ; iPaczvriski fc Sienkiewicz Il98ll ; 
iRappaport et al* 19831 ). The mass loss time-scale increases 
as M2 decreases, but the thermal time-scale of the secondary 
increases faster. When the thermal time-scale exceeds the 
mass-transfer time-scale, the secondary is not able to shrink 
fast enough (or in fact expands) in response to mass loss, so 
that the orbital evolution slowly moves back through longer 
periods. CVs in this phase of evolution are referred to as 
'period bouncers'. 

The rudiments of the theory of CV evolution, outlined 
above, have long been accepted, yet there are still signifi- 
cant discrepancies between the predictions of theory and the 
properties of the observed CV sample. One of the problems 
is that theoretical models persistently put the period mini- 
mum about 10 min short of the observed cut off (see iKolbl 
2002 for a review). Population synthesis studies also predict 
that ~ 99% of the intrinsic galactic CV population should 



consist of systems with P o r b S, 3 h, and that about 70% 
should be period bouncers (|Kolblll993l ; iHowell et ai]|l997l ). 
The first of these predictions is not reflected in the known 
sample of CVs — long-period systems comprise significantly 



more than 1% of observed CVs (see e.g. the left hand panel 
of Fig. [1} . Furthermore, on ly a handful of known systems are 
likely period bouncers (e.g JPatterson et al.ll2005l ). Clearly, if 
the standard theory of CV evolution is even approximately 
correct, the observed CV sample is not representative of the 
intrinsic population. However, since different types of CVs 
differ in intrinsic brightness, and all surveys are flux limited 
(and the majority severely so) there is no reason to think 
that an observed sample should be representative of the in- 
trinsic population. Apparent brightness and large amplitude 
variability are two factors that obviously affect the discov- 
ery probability of a CV. Since both intrinsic brightness and 
the frequency of large amplitude brightness variations de- 
crease with M, short-period, low-M systems are necessarily 
under-represented in the known sample. The question is only 
whether these and other selection effects can account for the 
size of the discrepancy between the period distributions of 
the predicted and observed samples. The right-hand panel 
of Fig. [T]goes some way towards illustrating the importance 
of the apparent magnitudes, by comparing the period distri- 
butions of bright and faint sub-samples of known CVs. The 
ratio of the number of short- to long-period CVs is larger 
when only apparantly faint systems are considered. 

Although detailed models of the intrinsic CV popula- 
tion exist, very little effort has been made to model even 
the simplest observational selection effects. It is therefore 
not clear to what extent theoretical predictions disagree with 
observations. Here we present a Monte Carlo approach that 
can be used to model any combination of observational se- 
lection effects quantitatively, and use it to demonstrate how 
common selection effects influence the observed CV period 
distribution. We then apply our method to a well-defined 
observational CV sample to show that the intrinsic CV pop- 
ulation predicted by the standard theory of CV evolution is 
indeed inconsistent with observations. 



2 THE COMPUTATIONAL METHOD 

We use a Monte Carlo model to determine the effect of ar- 
bitrary selection criteria on observationally constructed CV 
samples. The key parameters that describe a CV are Mi, 
M2, Porb and M. Given a description for the birthrate of 
CVs and a recipe for their subsequent evolution, population 
synthesis methods can be used to generate the predicted in- 
trinsic CV population as a function of these parameters. We 
generate a random sample of CVs drawn from such a dis- 
tribution, and distribute them within a model galaxy. By 
modelling the spectral energy distribution (SED) of each 
CV, we can construct a predicted sample corresponding to 
an observational CV sample subject to given selection crite- 
ria. 



2.1 The predicted intrinsic CV population 

The starting point of our modell ing is the p opulation of 
CVs resulting from model pm5 of lKolbl (|l993l ). This is the 
present day p opulatio n of C Vs that results from the birth 
rate model of Ide Kool (1 19921 ) when syste ms are evolved us- 
ing th e magnetic braking prescription of lVerbunt fc Zwaanl 
in the case of secondary stars with radiative cores, 
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and gravitational radiation alone for fully convective secon- 
daries The population synthesis model provides us with a 
probability distribution function (PDF) over the parameters 
Port, (M) (the secular average M), Mi, and M2. 

This model uses a bipolytropic treatment of the 
secondary, which becomes fully convective when M2 = 
0.25 Mq. The resulting period gap is between roughly 2.8 
and 3.2 h, which is not where the observed gap is. This 
shortcoming can in principle be overcome and does not have 
an important effect on the results we will derive, since the 
recalibration needed to improve the match to the observed 
location of the period gap has no serious repercussions for 
the overall population synthesis results. The PDF is zero at 
Porb = 2.92 h, and we use this value to distinguish between 
short- and long-period systems in the theoretical distribu- 
tion. Mi is kept constant throughout the evolution; i.e., all 
the material transferred from the secondary to the primary is 
assumed to be ejected during nova eruptions. Period bounce 
occurs at Mi ~ 0.07 MgQ, placing the period minimum at 
about 64 min. In this population, 99.7% of all systems are 
below the period gap, and 72% of all systems are period 
bouncers. 

Magnetic systems are not accounted for separately, so 
if they evolve differently, as has been proposed at least 
for polars (e.g. IWu 



irent l y, as nas been proposed at least 
& Wickramasinghe 19931; iLi et al. 1 11994 



IWebbink fc Wick ramasinghc 2002), they are not included 
in the population synthesis model. Furthermore, the largest 
secondary mass considered was 1 Mq , and all secondaries are 
assumed to be completely unevolved at the onset of mass 
transfer. 

It is worth pointing out that available population syn- 
thesis models assume a galactic age of 10 Gy or more. There 
is evidence that the galactic thin disc is only 7 to 8 Gy old, 
while the thick disc is no older than 10 Gy, with a scale 
height of ~750 pc, and a density mo re than a factor of 10 
lower than that of the thin d isc (e.g. iRevle fc RobinlfeoOll ; 
ISieeel et all |2002| ; lReidll2005h . If this is correct, the oldest 
CVs in the predicted population should be removed, which 
would reduce the predicted number of period bouncers. 



2.2 The model galaxy 

Our model of the Galaxy consists solely of an axisymmetric 
disk, with no spiral structure, no additional thick disc, no 
halo and no galactic bulge. This is consistent with the Popu- 
lation I assumption made in our population synthesis model 
and is reasonable, since the bulk of the intrin sic CV pop- 
ulatio n resides in the (thin) galactic disk (e.g. IStehle et al.l 
1 19971 ). 

Every system is assigned a position (r, z,0) in space, 
where z is height above the galactic plane, r is distance from 



1 iKolbl jl993h considered several combinations of birth rate mod- 
els and magnetic braking prescriptions, and found similar result- 
ing populations. The strength of magneti c braking i s, how ever, 
not known with great certainty, and llvanova fc Taaml j2003ft have 
pro posed a magnetic b r aking law that is significantly weaker than 
the lVerbunt fe Zw aan (1981) prescription. 

2 We will simply define the period bouncers in our population as 
systems with M2 $J 0.07 Mq. This is close to correct, although 
systems actually bounce at a small range of M2 ■ 



the galactic centre, and cj> is azimuth. The positions are gen- 
erated randomly in such a way that the overall distribution 
of CVs produces exponential density profiles: 

p(z) oc exp( — \z\/H z ), and (2) 

p(r) oc exp(~r/H r ), (3) 

where the radial scale length is H r = 3 000 pc, and the 
vertical scale height is 

{120pc for long-P or j, systems 
260 pc for normal short-P or f, systems (4) 
450 pc for period bouncers. 

Adopting different scale heights for these three sub- 
populations is motivated by the fact that the long- 
period, short-period (pre-period bounce), and period bounce 
regimes are successive stages in the evolution of a single 
CV. The age of a CV is the sum of the time spent as 
a detached binary (pre-CV) and the time since becoming 
semi-detached, and systems first make Roche-lobe contact 
at a range of orbital periods. This means that there is a 
spread in the ages of CVs at fixed period (besides that 
caused by period bounce), and, in fact, so me overlap be- 
tween the ages of sub-populations (see e.g. iKolb fc Stehlel 
1996). Nevertheless, sub-populations have different repre- 
sentative, average ages. Long period CVs are younger than 
~ 10 8 ' 4 y, while short-period, pre-period bounce systems are 
older than ~ 10 8 ' 8 y, and it takes ~ 10 9,4 y to reach the pe- 
riod minimum (e.g. iHowell et al]|200lD . The scale heights of 
the 3 sub-populations given in equation [4] were chosen to be 
representative of stellar populations with ages of ~ 10 8 ' 4 y, 
~ 10 9 ' 2 y, and ~ 10 9 6 y dRobin fc Crezelll986l ). 

We assume a ga lactic centre distance of 7620pc 
(Eisenhaucr et al. 2005h. and ignore the small offset of the 
Sun from the galactic plane; i.e., the Earth is placed at 
(r, z,(f>) = (7 620 pc, 0, 0). For computational simplicity, we 
divide the galaxy along the two planes of symmetry that 
include both the earth and the galactic centre (the planes 
z — and = 0) and simulate only one quadrant of the 
galaxy; i.e. z and ^ 4> < n. 

Interstellar extinction is computed by integrating the 
density of the interstellar medium 

7 620pc 



PISM K 6XP { 4 500pc ) X 6XP { 140^ ) (5) 

ijRobin et al.ll2003h along the line of sight to each system to 
obtain the neutral hydrogen column density, A^. We find 
the proportionality constant in equation [5] by setting Ay = 
30 for an object at the galactic centre. Extinctions in the 
U and .B -bands are taken as Au = 1.531Ay and As = 
1.324,4 v (|Rieke fc Lebofskvll 19851 ). 

By equation^ the G a ussian scale height of 190 pc found 
empirically by iPattersonl l|l984h is compatible with the ver- 
tical distribution we use for the youngest systems. This mea- 
surement of the scale height was based on the known sample, 
which is strongly biased towards intrinsically bright systems 
(this bias was even more pronounced two decades ago), and 
it seems reasonable that it may reflect the galactic distribu- 
tion of only the brighter, younger systems. Although younger 
stellar populations are expected to be more concentrated 
towards galactic mid-plane than older populations, the ob- 
serva tional evidence in the case of CVs is slim or non exis- 
tent (|Szkodv fc Howelll Il992l ; Ivan Paradiis et al.lll996l ). As 
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a check, we therefore also performed simulations using the 
same vertical density profile for all systems. In those simu- 
lations we used 

p(z) cc 0.80sech 2 (;z/323pc) + 0.20 exp(-|z|/656 pc), (6) 

which results from counts of M stars (|Gould et al.ll 19961 ). or 
a Gaussian profile, 

p(z)cxexp[-(\z\/H z ) 2 ] , (7) 
with scale height H z = 190 pc l|Pattersonlll984 l. 

2.3 The spectral energy distribution of CVs 

Four components, namely the accretion disc, bright spot, 
white dwarf, and secondary star, are included in the model 
of the overall SED for each CV in our simulations. 

We used a s light modification of the code described 
m iTvlendal (|l98lh to compute a 4D grid of accretion disc 
models over the variables Md (the rate at which mass flows 
through the disc), Mi, the orbital inclination i, and the outer 
disc radius Vd- We set the viscosity parameter a = 0.01 for 
M d «S 10 13 g/s, and a = 0.8 otherwis<0. The disc model 
assumes a steady state disc structure and produces a black- 
body spectrum for optically thick regions, while for a verti- 
cal optical dep th < 1, the same procedure as explained in 
ITvlendal ([1981) is adopted, so that emission lines appear in 
these regions. As will be explained in the next section, we 
take Md = (M) only for systems with stable discs (these are 
a small minority of the population). 

CVs are assumed to be randomly inclined with respect 
to Earth, so we select cos i for each system as a uniformly 
distributed random number between and 1. The inner disc 
radius is set equal to the radius of the primary (Ri), which 
is calculate d using the approx imate white dwarf mass- radius 
relation of iNauenberd (119721 ). The calculation of Vd will be 
explained in the next section. 

The bright spot is modelled as a blackbody extend- 
ing over an azimu th of 10° on the outer disc rim (see e.g. 
IWood et al.lll986h . 

The flux contribution of the second ary is computed 
using the spectral type -P or b relation of ISmith fc Dhillonl 
(|l998h : 

, . f 26.5 - 0.7P ort /h for P orb < 4h . 
P( ' \ 33.2 - 2.5P ort /h for P orb > 4h, 1 ' 

where Sp(2) = means a spectral type of GO, Sp(2) — 10 
is K0, Sp(2) = 20 is M0, etcQ. The secondary is assigned 
the absolute magnitudes that a main-sequence star of the 
same spectral type would have. This is not strictly correct, 
since CV secondaries are not in thermal equilibrium, but 

3 This choice was made mostly for computational convenience 
and is not very physical, but the influence o f a is small in all 
regards except U — B colours l lTvlendalll98l1) . This is discussed 
further in Section 4. 

4 Because the predicted period minimum is shorter than ob- 
served, we apply equation \E\ beyond the P orb range where it can 
be expected to be valid. Since these very low-mass secondaries 
are faint compared to other components of the binaries, this is 
not an important concern. The Mv max -P orb relation discussed 
below is also extrapolated to periods below the range where it is 
known to apply. 



unevolved secondaries are not expected ever to be large con- 
tributors to the optical flux. We may nevertheless improve 
our treatme nt of secondar y stars in the future by using the 
sequence of lKniggd l|2006h . In systems with M 2 < 0.07 M 
(i.e. period bouncers), the flux contribution of the secondary 
is neglected. 

We estimate the effective temperature (Te//) of the pri- 
mary using 




ijTownslev fc Bildstenl l2003h . The absolute U, B, and V 
magnitudes of the primary are the n found using the tab - 
ulation for pure H white dwarfs of iBergeron et all l|l995t ). 
In systems with Md > 2.0 x 10 17 g/s and i > 86°, we as- 
sume that the white dwarf is completely shielded by the disc, 
and makes no contribution to the flujJj. 



2.4 Dwarf nova outburst cycles 

Most CVs have values of M A that cause their discs to be 
subject to a thermal instability; these systems are the dwarf 
novae (we will abbreviate dwarf nova(e) to DN), and they 
spend most of their time with Md 7^ — M%, so that it is 
incorrect to use (M) in the disc luminosity calculation. DN 
are characterised by outbursts during which their luminosity 
increases by typically 2-5 mag. Outburst durations are days 
to tens of days, and the recurrence interval ranges from tens 
of days to decades. 

To determine whether a system is a DN, we use equa- 
tions 38 and 39 of lHameurv et ah! (|l998l ); systems with 

*>, 5 * 10 .v.» («)-(-,)"%, 

or 

have stable discs. None of the CVs in the theoretical intrinsic 
population satisfy the first inequality, so that we have no 
syst ems that are stable o n the lower branch of the S-curve 
(see lHameurv et al.lll998T ). All short-period systems in our 
model population are DN, and less than 0.1% of CVs in this 
population are permanently bright, nova-like (NL) systems. 

In order to calculate the disc spectrum correctly for the 
systems with unstable discs, we need to find the mass accre- 
tion rate through the disc in outburst (Mdo) and quiescence 
(MdCj)- These are related to the average value by 

(M) = CM d0 + (1 - C)M dQl 

where C is the outburst duty cycle (we assume that out- 
bursts are simple top-hat functions). For DN that have not 
yet reached the period minimum , we use the Mvmax-Porb 
relation of lHarrison et al.l l|2004h to predict the absolute 
V magnitude at maximum, but since the candidate period 

5 In detail, the issue of shielding by the disc is mo r e com - 
plicated than we have a ll owed for; see e.g. ISmakl il992l) , 
iMever fc Mever-HofmeisteJ dl982T) . and Knigge et al. (2000, 
2004). 
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Table 1. The DN outburst duty cycles we adopt. The final col- 
umn lists the fraction of the total intrinsic population represented 
by DN with each duty cycle. 





DN outburst duty 
duty cycle 


Percentage of the 
population 


Period bouncers 


5 x icr 3 


8.2 




5 x icr 4 


14 







20 


Normal DN 


0.1 


1.1 




0.01 


22 




5 x 10~ 3 


2.1 




5 x 10~ 4 


2.6 



bouncers reach brighter maxima than given by this relation 
(these systems show only superout bursts) , we use 



Mvmax = | 



5.92 - 0.383P or6 /h for normal DN 

5 for period bouncers. 



(9) 



We assume that all of the outburst light originates from 
the disc, and use a bo lometric correction of —1.8 (following 
IPatterson et al.ll2002l ) to find the disc luminosity (Ld) and 
hence the mass accretion rate through the disc at maximum, 



Mdo 



2R 1 Ld 
GMi 



(10) 



The rate at which mass flows through the disc in quiescence 
is then 



M dQ = 



(M) - CMao 
1 - C 



No clear empirical relation exists between the duty cycle 
and other observables, but for normal pre-period bounce 
DN a duty cycle of ~10% is in accord with observations 
(e.g. I Warner! [l995l ; lAk et all [2002^1. For period bouncers, a 
duty cycle of ~ 5 x 10 -3 is indicated by observations of the 
best candidate^. 

As it turns out, however, equation[5]and the mass trans- 
fer rates given by evolution theory cannot be reconciled 
with C = 0.1 (C = 5 x 10~ 3 ) for the majority of nor- 
mal DN (period bouncers). Duty cycles as large as these 
would imply (M) < CM d g and thus M d Q < 0. This is be- 
cause equation [10] gives Mdo between ~ 2 x 10 16 g/s and 
~8x 10 17 g/s, while (M) extends as low as 5 x 10 13 g/s. We 
deal with this problem by assuming that typical DN duty cy- 
cles are lower than those suggested by known systems. This 
is plausible, since the known sample of DN is likely to be 
biased towards systems with unusually frequent outbursts. 
Our grid of accretion disc models in fact extends only to 
Md = 3.2 x 10 12 g/s (by this point the disc is in any case 



6 Most (possibly all) short-period DN are SU UMa stars and 
have superoutbursts in addition to normal outbursts. However, 
since the latter are considerably more frequent, we ignore this 
complication. 

7 The duration and recurrence time of outbursts in these systems 
are uncertain. WZ Sge, EG Cnc, and AL Com are considered good 
perio d bouncer candidates (e.g. Patterson 1998; Patterson et al. 
2005), and appear to have outburst duty cycles of about 0.5%, 
but the observational baseline is not long enough for this to be 
firmly established. 



always fainter than the white dwarf). Therefore we require 

(M) 



M, 



dQ 



1 - C 



CMdo > 3.2 x 10 12 



For normal DN, if C = 0.1 does not satisfy the above in- 
equality, we try successively C = 0.01, C — 5 x 10 -3 , 
and C = 5 x 10~ 4 , and adopt the first value that gives 
MdQ ^ 3.2 x 10 12 g/s. For period bouncers, we start with 
C = 5 x 10~ 3 , and use C = 5 x 10" 4 if needed. If C = 5 x 10~ 4 
also leads to Md < 3.2 x 10 12 g/s, we assume that the sys- 
tem will never be observed in outburst, and that the disc 
luminosity is negligible in quiescence. Table [T] summarises 
the outburst duty cycles used. If we require only MdQ 0, 
then DN that must have C < 5 x 10 -4 comprise 19% of all 
systems, and less than 1% of the total population consists 
of DN that can have C > 0.1. 

Of course the procedure described above is very con- 
trived, but, given that it is not possible to reproduce the 
DN outburst properties expected on empirical grounds, this 
approach is at least reasonable. Most known DN were dis- 
covered through their outbursts. Assuming that (i) the space 
density of CVs is 5 x 10~ 5 pc~ 3 ; (ii) all DN reach M v = 5 in 
outburst; (iii) C = 5 x 10~ 4 implies an outburst recurrence 
time of 300 y (i.e. that 1/6 of these systems had an out- 
burst during the last 50 y); and (iv) every DN that reached 
my — 11 at maximum during the last 50 y was discov- 
ered and recognised as a DN, our treatment of outburst 
behaviour implies that we should know about 10 systems 
with C = 5 x 10 -4 . Given the uncertainty over space den- 
sity, the practically unquantifiable nature of the probability 
of discovering a DN in outburst, and the fact that we do 
know several DN that have had only one observed outburst 
(GW Lib being the best known example), we do not consider 
the — admittedly surprising — d uty cycles to be in themselves 
inconsistent with observations. IPattersonl (|l998l ) points out 
that a large population of low-M CVs could remain undis- 
covered if these systems have very long outburst recurrence 
times, and that, while this is not impossible, there is no 
empirical evidence of such behaviour. We can add that (ac- 
cepting the mass transfer rates that result from gravitational 
radiation alone, that DN outbursts are caused by an increase 
in Md, and that DN reach Mv — 5 at maximum) very long 
outburst recurrence times are not only possible, but indeed 
inevitable. 

We emphasise that equation is only used to estimate 
Mdo and MdQ, which are then used to calculate absolute 
magnitudes from our 4D grid of accretion disc models. We 
assume that the probability of catching a DN in outburst 
is equal to the duty cycle (which is correct for a single- 
epoch survey), and again use a Monte-Carlo method to sim- 
ulate some systems in outburst and others in quiescence. 
The outer disc radius is set to rd = 0.7Rl for NLs and 
DN in outburst, and Vd = 0.5Rl for quiescent systems 
|Harrop-Allin fc Warner||l996T l. where Rl is the Roche-lobe 
radius of the primary. 

Large sc atter in obser v ationally infe rred M at 
fixed period (IPattersonl 1 1984 IWarnerl Il987h . and — more 
compellingly — the coexistence of different subtypes of 
CVs in the same period intervals, imply that M also 
fluctuates around the secular mean on time-scales that 
are too long to be observed, but shorter than the bi- 
nary evolution time-scale (see e.g. IHameurv et alj Il989l ; 
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Table 2. Summary of the models considered here. 



Galactic structure 


DN outburst properties 


Model 


Equation [2] and [4] 


as described in 2.4 


Al 


Equation \E\ 


as described in 2.4 


A2 


Equation [7] 


as described in 2.4 


A3 


Equation [2] and [4] 


M d = (M) 


Bl 


Equation [6] 


M d = (M) 


B2 


Equation [7] 


M d = (M) 


B3 


Equation [2] and [4] 


M dQ =0 


CI 



Wu et al.|[l995l;lKing et al.ll 19951; [McCormick fc Frankll 19981 : 



Buning fc Ritterl |2004| ; iRitter et al.l l2000l ). The fact that 
these cycles are not included in the population synthesis 
model used here explains why our population contains no 
short-period NL variables. 

In addition to the treatment of DN outburst cycles de- 
scribed above, we will explore two limiting cases, namely 
(i) M d = (M) (which implies maximising M d Q and setting 
C = 0), and (ii) M d Q = (which is equivalent to adopting 
the highest allowed value of C). 

2.5 Summary of computational details 

In general, we find that the accretion disc is the brightest 
SED component in high-M CVs (i.e. DN in outburst and 
NL systems). The white dwarf dominates the optical flux 
in a significant fraction of DN in quiescence, both normal 
systems and period bouncers. The secondary star does not 
dominate the optical flux of any CV in our model. 

In order to allow for uncertainties in assigning appro- 
priate galactic scale heights and outburst duty cycles to our 
CVs, we have run a variety of models. All the models that 
we consider here are listed in Table [2] We view model Al 
to be the most realistic amongst these, and will focus on 
this model in the next section, but we also show results 
from other combinations of model galaxies and DN outburst 
properties for the sake of comparison (Section 4). 

We note again that our population synthesis model does 
not include evolved systems, and does not explicitly include 
magnetic CVs. Furthermore, our SED modelling does not 
make allowance for discless systems such as polars. These 
restrictions should be kept in mind when interpreting the 
results presented in the next section. 



3 SELECTION EFFECTS 

In the present study, we restrict our attention to three of 
the most common selection effects present in observational 
CV samples, namely magnitude limits, selection for blue 
U — B, and restricted ranges in galactic latitude. The im- 
pact of large amplitude variability (DN outbursts and nova 
eruptions), emission lines, and X-ray emission on discovery 
probability will be discussed elsewhere. 



3.1 Optical flux limits 

Bias towards finding apparently bright objects is the sim- 
plest selection effect, and a very important one. The proper- 
ties of (bolometric) magnitude limited samples have been 



studied before, with the assumptions that the observed 
radiation is the time-averaged accretion lumi nosity, and 
that CVs are dist ributed uniformly in space (Kolb] 1 19931 : 
iHowell et al. I 120011) . These assumptions make a compar- 
iso n with data almost meaningless (as was recognised 
by iKolbl 1 19931 ). Magnitude lim its together with galac - 
tic structure were con sider ed by Ritter fc Burkerj (|198(j ). 
IRitter fc Qzkanl l| 19861 ). and lRitterl l| 19861 ). However, the in- 
trinsic CV populations assumed in these studies were not 
based on full population synthesis models, only bright CVs 
(NLs and DN in outburst) were taken into account, and DN 
outbursts were not treated satisfactorily (it was assumed 
that all DN will be observed in outburst, and that Mi ex- 
ceeds M2 by a constant factor for all systems in outburst). 

Fig. [2] displays period histograms of samples produced 
by our standard model Al, with three different magnitude 
limits (V < 20, V < 16, and V < 14). In this figure (and, 
with the exception of Fig. [4] all other model period his- 
tograms we will show) the number of systems in each bin is 
scaled to reproduce a local space density of 5 x 10~ 5 pc -3 
(where local space density means counting systems inside 
a radius of 100 pejj. The shape of the period distribution 
depends strongly on the magnitude limit. As is expected 
from the fact that long-period CVs are intrinsically brighter, 
the fraction of the total sample made up by long-period 
systems increases for brighter magnitude limits; it is 4.7%, 
19%, and 33% for V < 20, 16, and 14 respectively. Dark 
grey shaded histograms represent period bouncers in each 
panel of Fig. [2] while the contributions of DN simulated in 
outburst are shown by the fine black histograms. 

Of course the known CV sample is not complete 
to 20th or even 16th magnitude, which raises the ques- 
tion of whether there is a magnitude limit for which 
we can compare known CVs directly to our predicted 
magnitude limited samples. Fig. [3] shows the apparent 
magnitude distr ibution of known non -magnetic, hydrogen- 
rich CVs from IRitter fc Kolbl |200o1) . with the distribu- 
tion in V of the complete model sample over-plotted in 
grey. Systems with P or t ^ 5h are omitted from both 
the observed and model populations, because CVs with 
evolved seco ndaries dominate the population at periods 
above ~5h jBeuermann et al.lll998l ; iBaraffe fc Kolbll2o"ool : 
iPodsiadlowski et al.ll200St ). The known samples with P or b < 
5 h and V < 12 and V < 13 contain only 5 and 17 systems re- 
spectively. A Kolmogorov-Smirnov (KS) test shows that the 
magnitude distribution of known CVs becomes marginally 
consistent with that of the model sample only at V < 13. 
The probability that the model and observed distributions 
are drawn from the same parent population is 0.08, 0.004, 
and 6 x 10~ 16 for limiting magnitudes of V < 13, V < 14 
and V < 16 respectively. This result is somewhat model- 
dependent. The inset in Fig. [3] shows cumulative distribu- 
tions for the observed (bold black) and model (dark grey) 



8 The C V space density is not known even to with i n a factor of 10 
fsee e.g.lRitter fc Burkertlll986tlShara, et aTlll986l: |Pohtandll996l ; 
iPattersod ll99St ISchwope et al] |2002t iThorstensen et al.l l2006t) . 

Our choice of 5 X 10 — 5 pc — 3 lies within the range suggested by 
different studies. When the model population is normalised to the 
formation rate of single white dwarfs, the predicted mid-plane CV 
space density is 1.8 X 10 -4 pc — 3 . 
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log(P orb /h) 

Figure 2. Magnitude limited samples produced by model Al (see Table [2) • The magnitude limit is given in each panel, together with 
the fraction of the sample with periods above the period gap. Here, and in all other predicted period distributions we will show, light grey 
shading of the range 2.63 h < P or t < 3.35 h indicates the model period gap (see Section 2.1). Bold histograms are all systems detected; 
the contribution of period bouncers is shaded in dark grey; the fine black histograms show DN that were simulated in outburst. The 
ratio of short- to long-period systems is 95:5, 81:19, and 67:33 at magnitude limits of V < 20, V < 16, and V < 14 respectively. Period 
distributions of all systems detected at each magnitude limit are plotted with linear vertical scales in the insets. 
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Figure 3. The distribution in apparent magni tude of known 
CVs (non-magnetic, hydrogen-rich systems in the iRitter fc Kolbl 
120031 sample; black), and a complete sample produced by model 
Al (grey; this histogram is arbitrarily scaled). Systems with 
Porb ^ 5h were omitted form both samples (see text). The inset 
shows cumulative distributions for known CVs (black) and the 
model sample (grey), together with a sample of objects that are 
uniformly distributed in space and all have the same absolute 
magnitude (fine black curve). Magn itudes used for the h istogram 
of the known sample are mag 1 from lRitter fc Kolbl l2003. and arc 
not in all cases V-band. 
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Figure 4. Orbital period histograms and cumulative histograms 
of the intrinsic population (dark grey) and samples with magni- 
tude limits of V = 20 (solid black) and V = 14 (dotted). The 
intrinsic and V < 20 distributions are both normalised, while 
the V < 14 histogram is scaled by the normalisation factor of 
the V < 20 distribution. The spike near the period minimum is 
increasingly suppressed at brighter magnitude limits. The broad 
feature at < log(P or (,/h) < 0.2 is strongly suppressed in the 
magnitude limited samples. 



samples. For comparison, the fine black curve represents ob- 
jects with identical intrinsic brightness, that are uniformly 
distributed in space. Since the known sample is likely to be 
incomplete even for V < 13, it should not be compared to 
the predicted magnitude limited samples presented here. 

The absence of a spike at the period minimum in the ob- 
served dist ribution has been the subj ect of several theo reti- 
cal papers. iBarker fc Kolbl (|2002h and lKing et~all (120021 ) aim 
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Figure 5. Period histograms of systems with V < 20 (dotted), 
and systems with V < 20 and the U — B indicated in each panel 
(black). To illustrate the effect of reddening, dark grey histograms 
show the distributions obtained when the same blue cuts are ap- 
plied to the intrinsic colours. Almost all systems with V < 20 also 
have both (U - B) < -0.7 and U - B < -0.7, and most long- 
period CVs have (U — B)o < —1.1. These samples were produced 
by model Al. 
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Figure 6. Period histograms of samples with the same U — B 
and (U — B)o cuts as in Fig. \E\ but here applied to a volume 
limited sample. Dotted histograms show all systems at d < 400 pc; 
black are samples with observed colours satisfying the blue cut 
indicated in each panel, and grey shows systems with intrinsic 
colours blue enough to satisfy the selection cuts. These samples 
were produced by model A3; i.e., the volume limited sample is 
representative of the intrinsic population. 



to explain this by removing the sp i ke fro m the intrinsic dis- 
tribution, whereas lKolb fc Bar affe (1999]) argue that it could 
possibly be a selection effect. Our results show that, relative 
to the whole population detected at a given magnitude limit, 
the spike at the period minimum indeed decreases in promi- 
nence at brighter magnitude limits, although it is present 
even in the sample with V < 14. Fig. [4] illustrates that the 
spike cannot be discerned in the intrinsic period distribu- 
tion. The broader feature between roughly log(P or ft/h) = 
and 0.2, caused by period bounce as well as the increasing 
evolutionary time scale of pre-period bounce systems, is very 
strongly suppressed, even at V < 20. 

3.2 Blue optical colours 

Most known CVs have (U-B) < -0.5 (e.g. lBruch fc Engell 
1 19941 '). This is because high-M systems have high continuum 
colour temperatures, while known low-M systems have the 
Balmer discontinuity in emission. However, in the systems 
with lowest M, the (accretion heated) white dwarf is ex- 
pected to dominate the optical flux, and this white dwarf 
may be quite cool (T e // < 10 000 K) . Selection cuts in U — B 
may therefore exclude the faintest CVs. 

Two UV-excess surveys covering large areas at high 
galactic latitude have produced reasonably large CV sam- 
ples. These are th e Palomar-Green ( PG) Survey (B < 16.2; 
U - B < -0.46; iGreen et al.lll986l ), and the Edinburgh- 
Cape Blue Object (EC) Survey (B < 16.5; U - B < -0.4; 
IStobie et allll997f ). Follow-up of the EC Survey is still far 
from complete; the PG Survey will be considered in more 
detail in Section [S] 

Fig. [5] shows the effect of a magnitude limit together 
with different U — B selection criteria. In every panel the 



dotted histogram shows all systems with V < 20, and the 
solid black histogram shows systems that also satisfy the 
blue cut. It is only at U — B < —0.9 that the blue selec- 
tion starts removing a large fraction of systems from the 
magnitude limited sample. 

In Fig. [6] we display P OT b distributions resulting from 
blue cuts imposed on a volume limited sample. It is seen 
that even U — B < —0.5 excludes systems near the period 
minimum, while U — B < —0.7 introduces a severe bias 
against short-period systems. This shows that many CVs in 
the intrinsic population are not bluer than U — B = —0.7, 
but, in the samples shown in Fig. [5j these systems were 
already excluded by the magnitude limit. Thus, UV-excess 
surveys that select objects for U — B < —0.7 are expected 
to be seriously biased against short-period, low-M CVs, but 
such a blue cut introduces hardly any additional bias in a 
survey that is also severely flux limited. 

The grey histograms in Fig. [5] and [6] display samples 
with the same blue cuts as shown by solid black histograms, 
but here the blue cuts are applied to the intrinsic colours to 
illustrate the influence of reddening. Notice that practically 
all long-period systems have (U — B)o < —0.9, but some 
get excluded from observed samples because they are on 
average more distant than short-period systems, and thus 
more highly reddened. This is amplified by the fact that 
they are (in model Al) also more concentrated towards the 
galactic plane. 



3.3 Restrictions in galactic latitude 

Surveys at high galactic latitude hold some promise of un- 
covering samples that are volume limited. In a galactic plane 
survey, on the other hand, intrinsically bright objects will be 
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Figure 7. The orbital period distributions of high galactic lati- 
tude (|6| > 30°; black) and galactic plane (|6| < 5°; dark grey) 
samples with magnitude limits V < 20 (left) and V < 16 (right). 
The plots are repeated for all our models of galactic structure. 
Surveys at different galactic latitude include very different popu- 
lations, and, as expected, the effect is larger for deeper samples. 
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Figure 8. Period and distance distributions of all objects in the 
volume defined by \b\ > 20° (grey) and systems in that volume 
with B < 20 (black). Only 40% of CVs in the volume limited 
sample are also in the magnitude limited sample, and practically 
all systems excluded by the magnitude limit have periods below 
the period gap, implying a strong bias against short-period sys- 
tems. Specifically, a B < 20 survey is incapable of finding most of 
the systems near the period minimum. This sample is from model 
A3, where the vertical density profile of stars is a Gaussian with 
a scale height of 190 pc. 

alistically small, the SDSS CV sample will by no means be 
volume limited. Clearly, even if no colour bias existed in the 
SDSS CV sample, modelling of selection effects is required 
to interpret it correctly. 

More fundamentally, accepting that CV sub- 
populations with different typical ages do not share 
the same vertical galactic distribution, a volume limited, 
high- 1 | sample will in any case not reflect the intrinsic 
galactic population. In fact, it is not practically possible to 
construct a volume limited observational sample in such a 
way that it represents the intrinsic population. 



detected out to larger distances than fainter systems. The 
direction of this effect is obvious; we simply show its mag- 
nitude in Fig. [7] with period histograms of high (black) and 
low (dark grey) galactic latitude samples with two magni- 
tude limits. Since this is a case where the results can be 
expected to depend sensitively on the assumptions we make 
about scale height, we include the corresponding histograms 
for model A2 and A3 (see Table [2} in Fig. [7] As expected, 
the ratio of long- to short-period systems is larger in the 
galactic plane samples in all cases; the effect is the smallest 
in model A2 (middle panels of Fig. [TJl , because this is the 
model in which the galaxy comes the closest to being an 
isotropic distribution of stars. 

In a survey covering the area |6| > 20°, only 4.1% of 
CVs in the survey volume are detected above a magnitude 
limit of B < 20 in the case of model Al. For model A3, 
40% of all systems in the survey volume have B < 20, and 
less than 0.5% of systems detected to this magnitude limit 
are long period systems; Fig. [8] illustrates this. The Sloan 
Digital Sky Survey (SDSS) reaches to g = 20, and has sky 
cove rage that can be roughly approximated as |6| > 20° 
(e.g. lYork etU]l2000h . Because CVs found in the SDSS will 
in future comprise the largest uniformly selected sample, 
this survey is bound to be very important in the context of 
obser vational studies of CV populations (e.g. ISzkodv et alJ 
2002). However, the result illustrated in Fig. [8] implies that, 
even assuming a galactic scale height that is probably unre- 



4 SENSITIVITY OF THE RESULTS TO OUR 
ASSUMPTIONS 

Since there are considerable uncertainties regarding the DN 
outburst properties of CVs and their distribution in the 
Galaxy, we explore here the effect of varying our assump- 
tions, by comparing several combinations of assumed galac- 
tic structure and outburst behaviour. In addition, we con- 
sider whether our results are sensitive to the viscosity pa- 
rameter we use in our disc SED models, since this is another 
potentially important source of error. 

At fainter magnitude limits, all our predicted magni- 
tude limited samples become more like the intrinsic distri- 
bution, and hence more alike. Therefore we consider here 
samples with bright limiting magnitudes. Fig. [9] displays pe- 
riod histograms of magnitude limited samples produced by 
the first six models listed in Table[2] The grey histograms are 
samples with V < 14, while black shows V < 16; the frac- 
tion of the total made up by long-period systems is given 
in each panel for both magnitude limits. This fraction does 
not depend very sensitively on the assumptions about galac- 
tic structure and DN outburst properties. As expected, it is 
more difficult to hide short-period systems if all CVs are 
assumed to accrete permanently at their long-term average 
rate, regardless of the adopted galaxy model. The fraction 
of long-period systems increases towards brighter magnitude 
limits in all cases. 



10 M.L. Pretorius, C. Knigge and U. Kolb 



Model CI produces the largest ratios of long- to short- 
period systems, and differs from, e.g., model Al also in pre- 
dicting fewer systems at a given magnitude limit, and in 
that relatively more detected systems are DN in outburst. 
Fig. 1101 compares the period distributions of all-sky samples 
with V < 16.5 resulting from model Al and CI. 

Because of the Balmer discontinuity, the U band is the 
region of the optical spectrum where the SED mo delling of 
the di sc emission is the least reliable. Figure 15 of iTvlendal 
l|l98ll ) shows that the U—B colour of low-Af discs is sensitive 
to the value of a (this sensitivity disappears at high M). For 
fixed inclination, varying a form 1 to 0.1 can change U — B 
by almost 0.35 mag for Md = 10 14 g/s and = 10 10 cm (the 
effect is larger for larger Vd, but the low-M discs in our model 
have small outer radii). In the very low-M systems where we 
take a — 0.01, the disc radiation has a negligible effect on 
the magnitudes and colours. Anticipating the application in 
Section we are for the moment only concerned with the 
influence that incorrect U — B colours may have when U — B 
cuts are applied to samples with bright magnitude limits. To 
determine whether an error of ~ 0.35 mag in the disc U — B 
in systems with 10 13 g/s < Md < 10 14 g/s is important, we 
performed the following test. The U — B colour of the disc is 
changed by adding 0.4 mag to the intrinsic disc U magnitude 
for all systems with M d < 2 x 10 14 g/s. Taking the B < 16.2 
sample produced by model Al, less than 0.5% of systems 
are lost from samples with blue cuts of U — B < — 0.7 and 
U — B < —0.9 when their original colours are altered in this 
way. The reasons are that low-M systems are rare in the 
magnitude limited sample, and that the flux contribution of 
the disc is increasingly unimportant compared to that of the 
white dwarf at lower accretion rates. 



5 APPLICATION TO THE PALOMAR-GREEN 
CV SAMPLE 

The CVs collected in the iRitter fc Kolbl (120031 ) catalogue 
include systems that were discovered on the basis of very 
different characteristics, such as physically distinct types of 
large amplitude variability, X-ray emission, blue colours, and 
emission lines. The resulting overall sample is too heteroge- 
neous to permit a reliable evaluation of its selection biases. 
Note also that, in addi tion to the discovery bias, a second 
bias is introduced in the lRitter fcPC olb (2003) sample by the 
need to find orbital periods. This is magnitude dependent, 
but it is also, e.g., easier to measure shorter periods (al- 
though this is not very important at P or b < 8 h) and there is 
a bias towards magnetic systems in deciding which CVs are 
worthy of follow-up observations to find P OT b- A quantitative 
comparison between theory and observations requires a rea- 
sonably large observed sample that has well defined selection 
criteria and quantified completeness. At the moment the ob- 
servational sample that conforms to these requirements the 
most closely is the sample of CVs found in the PC Survey. 

The PC Survey covered 10 700 deg 2 at [6| > 30° and 
was intended to have a magnitude limit of B < 16.16 and 
a col our selection criterion of U — B < —0.46 (|Green et alj 
1986). The blue cut and magnitude limit were always known 
to vary from field to field; the error in U — B and the com- 
pleteness were thought to be 0.39 mag and 84% respectively. 
More recent analyses of the PC data have shown that the 
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Figure 9. The orbital period distributions of magnitude limited 
samples produced by the first 6 models listed in Table [2] These 
are in every case for all-sky samples. Black (dark grey) histograms 
are for a magnitude limit of V = 16 (V = 14). In each panel, the 
contribution of long-period systems is given as a percentage of the 
total for the fainter (black) and brighter (dark grey) magnitude 
limits. 
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Figure 10. The orbital period distributions of samples produced 
by model Al (left) and CI (right). Both samples have a magnitude 
limit of V < 16.5. Dark grey shading and fine black histograms 
show the contributions of period bouncers and DN in outburst 
respectively. The contribution of long-period systems is given as 
a percentage of the total in each panel. 



actual magnitude limit and blue cut used by the survey were 
B < 16.16 and U - B < -0-71, wit h errors of a B = 0-34 
and a y-B = 0.24 (jjester et al.ll2005l : see also iLiebert et al] 
2005). The survey produced a sample of 36 CVs , of which 
3 are magnetic systems, and 6 have P or b > 5 h |Ringwaldl 
1 19931 ; Ringwald, private communication). Depending on the 
nature of V378 Peg (the only PC CV that still does not 
have a period measurement) the remaining 27 systems form 
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Figure 11. The orbital period distributions of the PG CV sample 
(including only non-magnetic system at P or b < 5 h; left hand 
panel), and the predicted sample (also with systems at P or b 
5h omitted; right hand panel). Light grey shading of the range 
2.2 h < P orb < 2.8 h in the left hand panel, and 2.63 h < P orb < 
3.35 h in the right hand panel indicate the observed and model 
period gaps. Dark grey shading represents period bouncers in the 
predicted sample. 

a statistically complete sample of unevolved, non-magnetic 

CVs. 

Following IJester et al.l (|2005l ). we take the probability 
of a system with B — B' and U — B — (U — B)' being 
detected in the PG Survey as P1P2, where 




We applied this detection efficiency to systems at |b| > 30° 
and P or b < 5 h in the population produced by our standard 
model Al. A period histogram of the resulting predicted 
PG CV population is shown in Fig. 1111 together with the 
distribution of the observed sample. 

It is evident from Fig. QT]that the distribution produced 
by the simulation is inconsistent with the observed one, the 
most obvious discrepancy being that the ratio of long- to 
short-period CVs is much higher in the observed sample 
than is predicted. If we conservatively assume that V378 
Peg is a short-period system, and taking the upper edge of 
the period gap as 3 h, 14 out of 27 observed systems (52%) 
at Porb < 5 h have periods above the gap. The prediction is 
only 6.7% (i.e. 2 systems). According to the binomial distri- 
bution, the probability that these numbers are consistent is 
less than 3 x 10~ 10 . We also note that none of the PG CVs 
is considered to be a candidate period bouncer, whereas the 
simulation predicts that 19% of the sample (i.e. 5 systems) 
should be period bouncers. 

It was shown in Section 4 that the properties of the 
predicted detected sample are fairly insensitive to poorly 
constrained aspects of our modelling of disc SEDs, galactic 
structure, and DN outburst properties. Of all the models 
we have considered, model CI allows for the largest popula- 
tion of short-period systems to remain undetected, because 
it combines a large galactic scale height for these systems 
with an extreme assumption of DN outburst behaviour. Re- 
peating this simulation using model CI yields a prediction of 
10% long-period systems in the PG sample, which is still in- 



consistent with the observed fraction. Based on this and also 
on the sheer size of the discrepancy between the predicted 
and observed PG samples, we are confident that the pre- 
dicted intrinsic population that is the basis of our modelling 
is ruled out by observations. We believe this is the strongest 
evidence yet that the standard disrupted magnetic braking 
model for CV evolution, with strong magnetic braking above 
the period gap and none below the gap, is incorrect. As has 
long been suspected, it predicts too many short-period sys- 
tems, and, specifically, too many period bouncers. 



6 CONCLUSIONS 

We have presented a Monte-Carlo technique that allows us 
to model selection effects in observed samples of CVs. Given 
any theoretical intrinsic population and any survey with 
well defined selection criteria, this technique can be applied 
to predict the observed population. We have explored the 
impact of flux limits, colour cuts, and observing restricted 
ranges of galactic latitude, and finally applied our method 
to the CV sample produced by the PG Survey. Our main 
results are the following. 

(i) The mass-transfer rates predicted by standard CV 
evolution theory are inconsistent with empirical estimates 
of the absolute magnitudes of DN in outburst, unless out- 
burst duty cycles are much lower than is generally assumed 
(i.e. < 10% for normal DN and < 0.5% for period bouncers). 

(ii) For practically achievable flux limits, no flux limited 
sample can be expected to reproduce the intrinsic CV pop- 
ulation. It is therefore quite inappropriate to compare an 
observed sample with the predicted intrinsic population. 

(iii) The properties of the period distribution of a mag- 
nitude limited sample depend strongly on the magnitude 
limit. The ratio of short- to long-period systems and the 
prominence of the spike near the period minimum decrease 
at brighter magnitude limits. 

(iv) The apparent magnitude distribution of CVs sug- 
gests that the known CV sample is not approximately 
magnitude-limited, even if a limiting magnitude as bright 
as V = 13 is adopted. 

(v) A blue selection cut of U — B < —0.7 introduces a 
serious bias against detecting short-period systems. How- 
ever, for existing UV-excess surveys, this bias is unimpor- 
tant compared to the effect of the survey flux limits. 

(vi) Magnitude limited surveys at high galactic latitude 
are expected to yield samples that are very different from 
those produced by galactic plane surveys, because they de- 
tect a larger fraction of all systems inside the volume defined 
by the galactic latitude range. However, a survey such as the 
SDSS is still not deep enough to be volume limited for CVs. 
This means that, for some time to come, the effect of a mag- 
nitude limit will have to be considered when comparing an 
observed sample to theory. 

(vii) A simulation of the selection effects present in the 
PG CV sample shows that observational biases are not suffi- 
cient to reconcile the intrinsic population predicted by stan- 
dard CV evolution theory (with strong magnetic braking 
above the period gap) with the observed PG sample. The 
real intrinsic CV population cannot contain as large a frac- 
tion of short-period systems (and, specifically, period bounc- 
ers) as is predicted by theory. 
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